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Mapping of the forces on biomolecules in cell membranes has spurred the development of effec¬ 
tive labels, e.g. organic fluorophores and nanoparticles, to track trajectories of single biomolecules. 
Standard methods use particular statistics, namely the mean square displacement, to analyze the 
underlying dynamics. Here, we introduce general inference methods to fully exploit information in 
the experimental trajectories, providing sharp estimates of the forces and the diffusion coefficients in 
membrane microdomains. Rapid and reliable convergence of the inference scheme is demonstrated 
on trajectories generated numerically. The method is then applied to infer forces and potentials act¬ 
ing on the receptor of the e-toxin labeled by lanthanide-ion nanoparticles. Our scheme is applicable 
to any labeled biomolecule and results show show its general relevance for membrane compartmen- 
t at ion. 

PACS numbers: 87.80.Nj, 02.50.Tt, 87.16.dp, 05.10.Gg 


The motion ofproteins and lipids in cell membranes and its relation to function have attracted considerable interest 
in recent years [l| . Motion is commonly followed by tracking of single biomolecules labeled by an organic fluorophore 
or an inorganic nanoparticle that allows detection via fluorescence, light scattering, etc. [Ij. Trajectories are usually 
analyzed by plotting the mean-square displacement (MSD) as a function of time. Parameters like diffusion coefficients 
and domain sizes are extracted by fitting MSD curves to analytical behaviors expected for different modes of motion, 
e.g. free Brownian diffusion, directed, confined or anomalous motion 0. 

A major physical motivation to biomolecule tracking stems from the actively debated origin of membrane com- 
mrtmentation. Free diffusion of membrane proteins in a sea of lipids was first postulated in the fluid mosaic model 
[sj. Following experimental observations of confinement, the lipid rafts Q and the picket and fence @ models were 
proposed. In the former, membrane proteins are preferentially located in domains with different lipid composition 
(lipid rafts). In the latter, compartmentation is ascribed to the combined action of the cytoskeleton and anchored 
transmembrane proteins, forming fences and pickets, respectively. Alternative models relying on more specific mech¬ 
anisms of protein-protein interactions have also been proposed [gI, 0 . Additional complexity arises from the fact that 
different confinement mechanisms may coexist and depend on the type of biomolecule [ 11 . 

The MSD-based approach has been used extensively. Alternative observables related to first-passage times @ or 
radial particle density distribution 0 have been proposed recently. More information on the dynamics is hidden in 
the full trajectory of biomolecules, though. Focusing on a single observable, e.g. the second-order moment for MSD, 
has the virtue of simplicity yet it wipes out information. In particular, it makes harder discriminating among different 
models of motion and does not provide systematic assessment of their validity. A more general approach based on 
inference methods [HI is taken by considering the likelihoods of the models themselves. A quantitative sense of their 
validity is thus obtained, together with systematic estimates of the parameters of the models and their uncertainties. 

Our aim here is to present a general inference approach to obtain maps of the forces and the potentials involved in 
the confined motion of biomolecules in cell membranes. Inferences are shown to provide sharp estimations of the local 
forces acting in the microdomains. We specifically consider the case of the receptors of e-toxins in the membrane of 
Madin-Darby canine kidney (MDCK) cells. 


The e-toxin is responsible for lethal enterotoxemia in livestock, due to the Gram-positive bacterium Clostridium 
perfringens (types B and D). A relatively inactive peptidic prototoxin is first synthesized and is then converted to a 
highly potent mature protein by cleavage and removal of terminal amino acids. The mature protein targets a specific 
receptor located preferentially in detergent-resistent domains of MDCK cells [l2|. The protein acts by heptamerizing, 
which leads to the formation of pores and the rapid modification of the membrane permeability to ions, causing cell 
death without any entry of the toxin into the cytosol in El. 

To label the e-toxin, we used 30-50 nm amine-coated lanthanide oxide nanoparticles (NPs) Yo. 6 Euo. 4 VO 4 (mean 
toxiniNP ratio, 1:1; see 0 ). These nanoparticles present several advantages: they are highly photostable without 
emission intermittency, they are synthesized directly in water and present extremely narrow emission, allowing efficient 
rejection of cell fluorescence 0 . Their size is directly determined from their luminosity 0 . Different emission colors 
are obtained using different lanthanide ions 0 . 
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FIG. 1: A 60-s trajectory of a Yo. 6 Euo. 4 VO 4 nanoparticle labeling e-prototoxin bound to its receptor on the membrane of an 
MDCK cell. The line color changes from red to blue (beginning/end of the trajectory). The motion is clearly confined. See 
EPAPS document No. for the movie at real speed (scale bar: 1 /xm). Excitation intensity, 0.2 kW/cm^; integration time, 50 
ms; readout time, 1.4 ms; r localization precision, 20 nm; temperature, 20°C. Inferred forces and potentials are shown in Fig. 
2A. The dashed lines indicate the mesh squares used for the inference. 


We used a wide-field inverted microscope (Zeiss Axiovert 100) equipped with a 63x, NA=1.4 oil-immersion objective 
and an EM-CCD (Roper Scientific QuantEM:512SC). The Eu^+ ions of the nanoparticles were excited with the 465.8- 
nm line of an Ar+-ion laser and their emission was detected using a 617/8M filter (Chroma). MDCK cells were grown 
to confiuency on glass coverslips. They were then rinsed, incubated with 0.04 nM of labeled e-toxin or prototoxin for 
20 min, rinsed 3 times, and observed in Hanks buffer containing 1% fetal calf serum and 1% penicillin-streptomycin 
either at 20 or 30°C. 

In all experiments (^400 cells), we observed several nanoparticles bound to a specific receptor on the cell membrane. 
We verified specificity of binding W pre-incubating the toxin for 1 h with an e-toxin antibody that prevents binding 
to the membrane (obtained as in [I 3 ) and verifying absence of nanoparticles bound to the cells. Toxins were kept at 
concentrations low enough to ensure that single toxins (and not oligomers) are tracked. Trajectories similar to those 
of prototoxins, which do not oligomerize, were indeed observed. The mean toximNP ratio 1:1 implies, assuming a 
Poisson distribution, that the fraction of NPs bound to zero, one and two or more toxins are 37%, 37% and 26%, 
respectively. Nanoparticles without toxins do not bind to the cells and are rinsed away. Given the size of the NPs, it 
is improbable that more than one toxin is present on the same area of the NP surface allowing simultaneous binding 
to more than one receptor. Furthermore, the binding ability of a fraction of the toxins may be impaired by the 
coupling to the NPs. We therefore estimate that the fraction of NPs bound to more than one receptor is less than 
10%. We also labeled e-toxins with the organic fluorophore Cy3 and observed again similar trajectories. This implies 
that the nanoparticle label does not modify the receptor motion, which is thus determined by the receptor mass and 
the membrane characteristics (viscosity, forces, etc.). The receptor motion was studied during 150 to 300 s. Figure 1 
shows a portion of the confined trajectory of a prototoxin bound to its receptor. We verified that, given the diffusion 
coefficient and the domain size, we are not limited by the image acquisition time (21.4 or 51.4ms) [l9|, l20|. Relatively 
short portions of trajectories were considered, so as to exclude possible drifting of the membrane domain, the cell or 
the microscope setup. 

Langevin equations for the position r{t) and the velocity v{t) of a biomolecule subject to molecular diffusion and 
to the force induced by a potential V are: 


dr 

dt 


dv 


= V ; (^) + . 


( 1 ) 


Here, m is the mass of the biomolecule, 7 and D are the friction and the diffusion coefficients inside the microdomain 
of the membrane. The zero-average Gaussian noise rapidly fluctuates in time, accounting for the effect of thermal 
noise. Smoluchowski’s overdamped approximation [ 2 l| to ([T]) is sufficient for the motion of biomolecules. Indeed, 


the typical time for the relaxation of the velocity to local equilibrium isr= — ^10 (since m 10 kg and 
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FIG. 2: Force and potential fields inferred inside two membrane microdomains. The length of the arrows is proportional to the 
magnitude of the force. The potential is plotted as level curves with the bar indicating the potential’s amplitude on the isoline. 
The interpolation scheme described in the text was used, leading to an effective time step of 5 ms. The diffusion coefficent 
for A is 4.75 x fim?s~^ and for B is 8.15 x . Posterior distributions for the five locations indicated by the 

squares are shown in Fig. (4] 


7 10 ^ kg/s [HI). Hence, the velocity is slaved to its local forcing and ([T]) reduces to 

at 7 

The Fokker-Planck equation [ 2 ^ associated to © reads 


( 2 ) 


dtP =--V ■ (F P) + DAP, (3) 

7 

where the force F = — VH. Kolmogorov equation © governs the transition probability P(r,t|ro,to) to get to the 
space-time point (r’,t) conditional to the initial space-time position (ro,to) of the bio molecule. It follows from ([3]) 
that the probability P can be expressed as a path integral [2^ over all paths r{s) connecting ro to r: 


P(r, t|ro, to) oc J P)r{s)e (4) 

The term Q{r{s)) = {dr{s)/ds — /AD is the quadratic Gaussian weight governing the probability of 

displacements over an infinitesimal time interval. 
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FIG. 3: Typical evolution of the inferred value of a local force (left) and its standard deviation (right) with the number of 
points used to infer them. Note the rapid convergence to the real value of the force (indicated by the solid horizontal line). 


In practice, space is discretized in a fine regular mesh of n? squares, as shown in Fig. 1 (n = 8). The size of the 
mesh is taken small enough for the forces to be smooth on that scale. At the lowest order, forces are approximated by 
a constant value within each mesh square and we shall show later that higher-order variations are indeed negligible. 
The integral appearing at the exponential in (|1]) is approximated by the corresponding discrete Riemann sum (see 
dni). Mesh squares are indexed by the pair (i,jf) (with i,jf = 1,... ,n) and the force acting in is denoted 
Fi^j. Our goal is to estimate the 2n^ + 1 unknowns U = {D^{Fi^j}}^ i.e. the forces and the diffusivity within a 
membrane subdomain, governing the trajectories of the labeled biomolecule. Variations in D can be handled similarly 
(see below). 

Inferences methods (see, e.g., [HI) generally feature two steps: a) the derivation of the posterior probability 
distribution of the unknown parameters of the model given the experimental observations; b) sampling from the 
posterior distribution to estimate the parameters. Specifically, it follows from Bayes rule that the posterior probability 
distribution P {U\T) of the set of unknown parameters U given an observed trajectory T reads 

, 5 ) 


where P {T\U) is the likelihood of a trajectory given the parameters U and P (T) is a normalizing constant. Po{U) is 
the prior probability, which we take constant. As for the sampling part (b), we used Monte Carlo methods to compute 
the average over the posterior distributions. The latter are generally well-peaked and maximum values provide then 
good estimates of the average values. 

An asset of our specific problem is that the diffusivity D is the only global parameter whilst the forces Fij 
appear in the likelihood additively at the exponential. It follows that the contributions of the various squares of the 
mesh factorize as P {U\T) = P {Fij, D\T). The contribution of each mesh square reads 


exp 

P{Fij,D\T)<x n - 


4DAt 


AirDAt 


(6) 


Here, /i indexes the various time steps (discretized by At) and the product is restricted to those times when the 
biomolecule is detected within the mesh square Sij. Note that discretization introduces a fortiori an ambiguity when 
the bio molecule crosses the lines of the mesh and moves to a new square. The choice made in is to simply use 
indices of the starting square. Corrections will be shown shortly to be negligible. 

The crucial element ensuring well-peaked posterior distributions and sharp inferences is that the trajectories of the 
biomolecules are well confined to subdomains. It follows that the various squares of the mesh are crossed multiple 
times, permitting the acquisition of a massive amount of information. Even for those squares where the largest forces 
are measured, i.e. the residence time is the shortest, the amount of data is sufficient to permit sharp inferences. Note 
also that posterior distributions for the forces are Gaussian, as seen directly in ®. 

To have a quantitative sense of the quality of the inference scheme, we numerically generated ensembles of trajecto¬ 
ries with the same force fields and diffusion coefficients as obtained from the experimental data. Posterior distributions 
were found to be sharply peaked at the values used to generate the trajectories. Typical evolutions of inferred values 
vs the number of points used for the inference are shown in Fig. [3l Convergence is manifestly rapid and the standard 
deviation brackets the real value even for few data points, providing a sensible estimate of the error bars. Predictions 
by our method were found to be more precise and to require less data points than those based on a single statistic. 
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FIG. 4: Posterior probability distributions of the forces at the locations indicated with the squares in Fig. O The curves for 
horizontal (vertical) components of the force are 2,4,5 (1,3). Solid curves are obtained by the interpolation scheme described in 
the text. Note that average values are extremely close to those obtained without interpolation (dashed curves) yet the variance 
is reduced. 

e.g. MSD or radial particle density distribution. In summary, simulations provide strong support to the validity of 
the inference method. 

To visualize the results, it is convenient to plot the potentials V as in Fig. [2l To that purpose, the potential is 
written as a polynomial of order C (C = 4 in Fig. [2j) : V (r) = Yl'j=oYli=oCtijry^ \ The constants are fitted 
to the experimental force fields, minimizing the squared error by standard simplex methods. Potentials that we find 
are incompatible with a cytoskeleton fence-type model, where a steep wall-like potential is expected. This type of 
domain, however, may still influence the receptor trajectories on time scales below our resolution. Variations within 
microdomains for the diffusivities were found to be small, i.e. about 4% vs 65% for the forces. 

Discretization errors were controlled by the following method. Given two acquisitions (xi, G) and (X 2 G 2 = G + At), 
we interpolate their transition probability by summing over all possible positions x' at the intermediate time t' = 
ti-h At/2, i.e. P (^ 2 , t 2 |xi, ti) = / dx'P (x 2 , t 2 |x', t') P (x', t'l^i, ti). The process can be further refined by introducing 
additional intermediate points. The effect of the interpolation mostly amounts to a reduction of the error bars, without 
any major shift in the estimates of the forces, as can be seen in Fig. IH 

In conclusion, we have developed an inference approach that fully exploits information hidden in labeled biomolecule 
trajectories. The technique is generally applicable to any type of biomolecule and trajectory, including intermittent 
trajectories like those obtained with blinking quantum dots, and for forces and domains that change in time. We have 
explicitly demonstrated the value of the method by mapping the forces and the potentials involved in the confined 
motion of the e-toxin receptor in the membrane of MDCK cells. Results obtained here indicate that the method, 
especially in combination with data on cytoskeleton destruction and cholesterol depletion, is poised to shed light onto 
the controversial mechanisms of membrane compartmentation. 
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